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Abstract 

We define the harmonic evolution of states of a graph by iterative appli- 
cation of the harmonic operator (Laplacian over Z2). This provides graphs 
with a new geometric context and leads to a new tool to analyze them. The 
digraphs of evolutions are analyzed and classified. This construction can 
also be viewed as a certain topological generalization of cellular automata. 

MSG 05C50, 15A33, 05C75, 05C85. 

1 Introduction 

The Laplacian of a graph is a well known tool to study graphs. Typically, we 
retrieve its eigenvalues and deduce certain graph's properties [1, 2, 5]. In this 
paper we rather analyze the graph's "harmonic operator", the Laplacian mod- 
ulo Z 2 , or — more concretely — the monoid it generates. Interestingly, the 
emerging algebraic properties may be interpreted pictorially and visualized. 
We define a cellular automaton on an arbitrary (simple) graph as follows. 

The Harmonic Game ("game of light"). Choose a graph to represent your 
"universe". Each vertex may be in one of the two states: "on" or "off" (excited 
or not excited). Any initial state of a graph will evolve stepwise according to 
certain "laws": 

1. If a non-excited cell has an odd number of excited neighbors, it becomes 
excited. 

2. If an excited cell has an odd number of non-excited neighbors, it will 
remain excited, otherwise it relaxes and goes to rest. 

Examples of evolutions on various graphs may be inspected in Figures 1 and 2. 
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Figure 1: Examples of harmonic evolution on graphs ("game of light"). 

Clearly, each process has to eventually end up in a loop, consisting possibly 
of a single element, like Fig. lb. Even small change in structure may radically 
change the size of the loop and of the path towards the loop (cf. Fig. la and lb). 

The question is — what are all possible evolutions on a given graph G due 
to different initial states. Laborious exercise would recover a digraph whose 
nodes are the states (2™ of them, n=number of vertices in G), and whose di- 
rected edges indicate the consecutive states. See Figure 3 for an example. It 
will be called the evolution digraph of G and denoted G*. It turns out that 
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Figure 2: Game is like a propagation (a) a string with 9 nodes (b) a diamond 
shaped graph. 

evolution digraphs have the same general structure: each G* consists of a cer- 
tain number of closed loops, and to every node in a loop a certain tree is at- 
tached, constituted by states descending towards it. Interestingly, every tree in 
G* has the same form. Other examples of graphs and their evolution digraphs 
are shown in Figure 5. 



The structure of these trees and of the loops is concealed in the harmonic oper- 
ator and its powers. We will describe an algorithm that uncovers it. 

The harmonic game reminds one of Conway's "game of life." But things are 
here more fundamental. Unlike the usual cellular automata, the rules are not 
arbitrary, but are derived from the topologically invariant notions, that of bound - 
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Figure 3: Evolution digraph G* of a graph G. 



ary and coboundary. They represent the action of the harmonic operator a. 
Moreover, the rules and the processes in Figures 1 and 2 present only the sur- 
face of the actual dynamic, because both the vertices and the edges are actually 
involved in the process. 

Motivations for studying Harmonic games (HG): 

1. (Graph theory) A tool to analyze graphs in an alternative way to that of 
Laplacian's eigenvalues. The behavior of the states convey information on 
eigenspaces of powers of a. In a sense the harmonic game is a method to il- 
lustrate the structure of eigenspaces of the Laplacian by trees, loops, etc. 

2. (CA) Although we see it as a graph-theoretic study, we view HG as a step 
beyond the regularity of grids used in cellular automata. And yet another dif- 
ference: HG conceals actually two intertwined evolutions, one on vertices and 
one on edges, as explained in the next section (the rules presented here present 
only a surface of the dynamics). 

3. (Metaphor) A toy model of reality. HG is a metaphor of a typical theory 
of theoretical physics. (The graph itself is like a protogeometry: locations and 
connections only). One is reminded of Wheeler's search for the "Law without 
Law" [6] or his "pre-geometry" [7]. 

Remark: If biologists have their mathematical toy in the "game of life", HG 
may be viewed as a physicist's toy, a metaphor for "light propagation": The 
rules for a one-step change defined by the harmonic operator are styled after 
the harmonic operator in differential geometry (and Maxwell's theory of light). 
Harmonic evolution has indeed the flavor of light "propagation" (see Fig. 2). 
It actually obeys "Huygens principle": the sum (mod 2) of two evolutions is 
an evolution. Note other features: the superposition principle, periodicity in- 
duced by boundary conditions, duality in propagation (like in electromagnetic 
wave), simplicity of laws, their topological foundations, etc. 
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2 Harmonic Evolution 



Let G = (Ver (G), Edg (G)) be a simple graph (no loops, no multiple edges [5]) 
with the set of vertices Ver(G) and the set of edges Edg(G). We define a state 
(respectively co-state ) as a subset of vertices (respectively edges) of G. Denote 
the family of all subsets of Ver(G) by V G , and that of Edg(G) and by E G . Define 
addition as the symmetric difference A + B =: (A — B)\J (B — A) for any two 
subsets A and B ("adding modulo two"). 

The boundary operator 8 : V G — > E G sends maps into co-states, namely for any 
S G V G , we define 8S as the sum of the sets of edges adjacent to vertices of S, 
added modulo two. Similarly, a co-boundary operator 8 : E G — > V G maps co- 
states into states, namely for a G E G we define 8a as the set of vertices adjacent 
to the edges of a, added modulo two (see Figure 4). 




Figure 4: Boundary 8, co-boundary 8, and harmonic map a = 88 

By analogy to Hodge theory in differential geometry, the harmonic map is de- 
fined as the composition a = 88 + 88. In particular, if restricted to states, the 
harmonic operator becomes 

a = 88 e EndV G (2.1) 

By iterative application of the harmonic operator to an initial state S G V G , we 
obtain a sequence {S, aS, a 2 S, . . .} that we call the evolution of S, examples of 
which are in Fig. 1 and 2. 

For a given graph G, all evolutions form an evolution digraph G* , the nodes 
of which are the states, V G , and the oriented edges of which are of the form 
(S, aS) (see Fig. 3). The problem is to determine the evolution digraph of a 
given graph. 

It turns out that digraph G* consists of certain number of cycles, (later collec- 
tively called the loop ensemble), to each node of which attached is a copy of a 
descending tree (later called the characteristic tree of G). Digraph G* is a source 
of invariants of G. 

Before we go on, inspect Figure 3 and 5 that present evolution digraphs for a 
few small graphs (ignore for now the labels in Figure 5). Notice these features: 
(i) the sum of nodes in the loops is a power of two; (ii) the size of each loop 
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Figure 5: Three graphs and their evolution digraphs. 

divides the size of the longest loop; (iii) each descending tree has the same 
structure; (iv) the number of nodes in each tree is a power of two; (v) trees have 
some regularity in shape. Thus the shape of the tree and of the loop part would 
suffice to reconstruct the digraph, as the arc orientations are unambiguous. 

3 Geometry of Graphs 

The following 'geometrization of graphs' will be exploited. The space of states 
Vg of a graph will be viewed as an n-dimensional linear space over Z 2 , that is, 
Vg = Z£, where n = |Ver(G) |. The structure of the graph is represented by two 
Z-valued scalar products in the space Vg. One is the natural scalar product g, 
which for the vectors representing single vertices is g(vi, vj) = Sij. The other is 
a (possibly degenerated) scalar product A represented in the natural basis by 
the adjacency matrix, that is A(vi, Vj) = 1 if (v u Vj) e Ver(G), and A(vi, Vj) = 
otherwise. 

The harmonic operator (2.1) is geometrically an endomorphism of space Vg, 
self-adjoint with respect to g. One may easily show that in the natural basis of 
vertices, it has the matrix form 

a = A+D , 

where D is a diagonal matrix with values 1 at the diagonal entries that corre- 
spond to the odd vertices, and otherwise. This is the Laplacian of the graph 
taken modulo 2. As a matrix, a has an even number of l's in every row and 
every column. Let fl denote a column of all entries equal to one. Then actu- 



5 



ally any square matrix with entries and 1 only satisfying (i) a T = a and (ii) 
aO = 0, determines a graph. 

Evolution of a state S e V is given by S(t) = a l S where (eN plays the role 
of time. Thus, in order to study harmonic evolutions we need to start with 
the monoid of endomorphisms generated by a. Due to the finiteness of G, the 
monoid is of the form: 

M(a) = {1, a, a 2 , . . . , a k , . . . , a k+n ^} , 

where a k+n = a k for some k and n. Part L(a) —: {a k , . . . , a k+n ^ 1 } forms a 
cyclic group, part T(a) =: {1, . . . , a k } is called the tail of M(a). We shall use the 
notation n = \L(a) \ and k = \T(a)\. 

Lemma 3.1 Operator tt = a v with p — \L(a) \ ■ \T(a) \ is an orthogonal projection in 

V. 

Proof: One may easily show that n is the neutral element of the group L(a). 
Therefore, 7r 2 = n, i.e., ir is a projection. Orthogonality of ir follows from self- 
adjointness of a. □ 



Proposition 3.2 The space of states V G is a direct product of two subspaces 

V G = T(G) 8 L(G) (3.1) 

where T(G) — Kern and L{G) = Ini7r. Moreover, restriction of a to T(G) is nilpo- 
tent, a k = 0, and restriction of a to L(G) is an automorphism ofY. 

Lemma 3.1 implies: 

Corollary 3.3 The dimension dim L(G) is even for any graph G. 

Proof. Indeed, projection n must have two eigenvalues: namely 1 for the 
subspace L(G), and for T{G). Therefore, a basis of Vg exists in which tt is 
expressed by a diagonal matrix tt = diagfl, 1, . . . , 0, . . . , 0]. Trace does not 
depend on basis; therefore, dim L = Tr 7r = Tr it. Since Tr a 1 is even for any i, 
so is dim L. rj 

Proposition 3.2 explains the general structure of evolution digraphs. Here is its 
meaning. 

(z) The n-cube Vg = may be considered as a digraph (Vg, c H), where the 
"looped arrow" denotes the relation of succession. That is, for any two vectors 
(states) x,y e V G we say x °r> y only if y = ax. 

(ii) Subspace T(G) € Vg is closed under (nilpotent) action of a. It may be, too, 
considered as a digraph (T(G), S->), with a similarly defined relation of succes- 
sion. Due to nilpotency of a, this digraph has a form of a tree with as the root. 
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We shall call this digraph the characteristic tree of G. 

(Hi) Similarly, the subspace L(G) e Vg is closed under action of a, but this 
time a acts as an automorphism. Thus, as a digraph (L(G),<++), it consists of a 
number of loops (cycles) arising as the orbits of the action of the cyclic group 
generated by a. We shall call this digraph the loop ensemble of G. 

(iv) Now, the evolution digraph G* results as a "semidirect" product of the two, 

G* = L(G) x T(G) = {L(G) x T(G),<^} 

with the succession inLxT defined 

(x, y) °^ (x 1 , y') if either (x = x' and y R-> y') or (x R-> x' and y = y' = 0) . 

As a digraph, G* = T(G) x L(G) will be simply denoted as G* = {T(G),L(G)}. 

In the following two sections we shall characterize the two digraphs. 



4 Classification of Harmonic Trees. 

Let A * B denote the digraph obtained as the direct product of two digraphs A 
and B. That is we set (a, b) -< (a', b') iff a -< b in A and b -< b' in B. 

Theorem 4.1 The characteristic tree of a graph can he uniquely factored into a product 
of binomial trees: 

T(G) = il 1 * 4 2 * . . . * l\ k (4.1) 
where Ii denotes a binomial tree of height i. 

Proof: Since a restricted to T(G) is nilpotent, then — by Jordan decomposi- 
tion theorem — there exists a basis in T(G) such that operator a takes in it a 
quasi-diagonal form: 



Ji 



J,. 



where J, = 





1 
1 







(4.2) 



(possibly, some J, are lxl null matrices). This means that the space T(G) 
decomposes into eigenspaces of the harmonic operator, on which map a acts 
as a "raising" operator. Namely, denote the vectors that span the subspace E{ 
corresponding to the sub-matrix Jj by {/i, . . . , f q }. Than a acts : 



fi 



fi 



and 



(4.3) 
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Figure 6: Products of binary trees: a fragment of multiplication table. 

Note that each /, (for i > 1) is an image of exactly two states, namely and 
fq + fi-i- So is, therefore, any sum of these vectors. The state fx, and any sum 
J2iei fi containing fx, do not have any anti -images through a. Therefore, the 
states of the subspace Ei form a binomial tree of height q with 2 q edges, among 
which half, 2 9_1 , form the top (starting) nodes. Decomposition of a into lin- 
early independent blocks (4.2) corresponds to multiplication (4.1) of such trees. 

□ 

A few products of binary trees are shown in Fig. 6. Note that the height of the 
resulting tree is that of the highest tree in the product. 

The above theorem implies that the number of the binomial trees in the de- 
composition (4.1) equals dimKcra. Indeed, in the quasi-diagonal form (4.2) of 
the harmonic matrix, each Ji contributes a single one-dimensional subspace to 
Kera. (Clearly dimKera \t(g) = dimKera, since Kera C T(G)). Another im- 
plication is that if dimKera = 1, then T(G) is a binomial tree, T(G) = I\T(a)\- 
If dimKera = 2, then T(G) = I\t{o)\ * IdimT{G)-\T{a)\- 

The composition of the characteristic tree, that is the exponents of the factor- 
ization (4.1), may be deduced by means of the elements of the monoid M(a). 

Theorem 4.2 The multiplicity of the binomial tree Ij in T(G) is given by the follow- 
ing formula 

rij = 2dim Kera^ ' — dim Kera J+1 — dim Kera 1 ^ for j < k (4.4) 

Proof: The dimensions of kernel of powers of a satisfy the following system 
of equations: 



dim Ker a 
dim Ker a 2 
dim Ker a 3 



= rii + n 2 + ■ ■ ■ + n k 

= rix + 2(n 2 + . . . + n k ) 

= nx + 2n 2 + 3(n 3 + n 4 + . . . + n k ) 



(4.5) 



dim Ker a k 



ni + 2n 2 + 3n 3 + . . . + kn k 
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Denote Ki = dim Ker a*. Then (4.5) is equivalent to a matrix equation K = Mn, 
where K'Uj = min(i, j). The inverse matrix M _1 , that is easy to find (here for 
k=5 for illustration): 



11111 

1 2 2 2 2 

1 2 3 3 3 

1 2 3 4 4 

1 2 3 4 5 



— 1 



r 2 -1 OO-i 

-1 2-1 



K = 



K 



0-1 2-10 



0-1 2 
L -1 2 J 



solves (4.5) into (4.4), as claimed. ^ 

This theorem gives the algorithm for retrieving the structure of the character- 
istic tree from the harmonic matrix a of a given graph. 

5 Harmonic Loops 

Now we shall look for the structure of the loop ensemble of a given graph. 
The elements of group L(a) act — if restricted to the subspace L(G) e Vg — 
as automorphisms. The action of L(a) decomposes the subspace L(G) into a 
number of orbits. Since the group L(a) is cyclic, the orbits are loops. (Clearly, 
the orbits do not form linear subspaces of L(G)). It follows thus: 

Proposition 5.1 The least common multiple of the lengths of the loops in L(G) is 
equal to the order of the loop group, m = \L(a)\. There is a loop in L{G) of length m. 
In particular, the rank of any vector of L( G) divides m. 

Denote a loop of length i by Lj. If nLi denotes a disjoint sum of n loops of order 
i, then the decomposition of L(G) into loops may be expressed symbolically as 
a formal sum 



(for examples this notation see Fig. 5) Due to Proposition 5.1, the only nonzero 
numbers could be those with i dividing m = \L(a)\. The total number of 
loops is n = J2 n f Clearly all lengths of the loops sum up to a power of two, 
namely £. i- n, = \L(G) \ = 2 dimi ( G ). 

Now we recover the exact structure of the loop ensemble. That is, we retrieve 
the coefficients of the decomposition (5.1), from the monoid M(a). We shall 
denote the eigenspace of an endomorphism g by F(g), and its dimension by 

<lim /' I'/:. 

Theorem 5.2 The number of loops of length p in L(G) is 



L(G) = n\L\ + n 2 L 2 + . . . + n m L. 



'in 



(5.1) 




(5.2) 



where fj, denotes Moebius function, and the sum runs over divisors of p. 
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Proof: First, notice that F(a % ) C F(a j ) only if i\j. Hence, for any j we have 

3 3 

\F(a j )\ =Y,d-n d = Y,d-n d -ii{d,j) (5.3) 



d=l 



where ^ is the "factorization function" defined: /J-(i,j) = 1 if and /x(i, j) = 
otherwise. For conciseness, denoted = |i^(a*)| = 2 dimF ( ai ) (this is the number 
of states in the eigenspace of a 1 ). The system of equations (5.3) has form F = 
An, where A is a lower-triangular matrix 



1 

1 2 

1 3 

12 4 

1 5 

1 2 3 6 



or 



Aij = i ■ fi{i,j) 



(5.4) 



It is well-known that the inverse of matrix M*j = fi(i,j) is given by (M 1 ) i 
fj,(j), where \i is the Moebius function, defined 



M (d) 



(-ir 





if d= 1 

if d = Pl • p 2 ' • • • ' Pr 

if p 2 \d for any prime p 



where p\ , p 2 , ■ ■ ■ , p r are distinct primes (for the properties of the Moebius func- 
tion see [4]). Equation (5.4) is a modified version of the factorization function, 
and it easily solves to (5.2). q 

The above theorem provides an algorithm to construct the structure of the loop 
ensemble of the evolution digraph of a given graph. 

Remark 5.3 Here is alternative view on the inversion (5.2) of the system (5.3). 
Let {N, ^} be the poset (lattice) of natural numbers with the order induced by 
multiplicative relations, i.e. a -< b if a \ b. Let [a] denote the sub-poset of all 
elements below a, i.e. [a] = {n e N | n -< a} (a filter generated by a). Let a be 
the set of maximal elements in [a] — {a}; say a = {a\, . . . , a n }. The subsets [a{\, 
[a 2 ],...,[a„] form a "H-algebra", i.e. their mutual intersections are again some 
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Figure 7: Examples of graphs and their loop and tree structure. 



subposets generated by elements of 



30 



1 



60 



/ \ 30 20 12 

10 15 / I X I X I 

| \ / | 15 10 6 4 

2 5 3 1X1X1/ 



5 3 2 

\ I / 
1 



The numbers of vertices in loops of order k may be obtained by the inclusion- 
exclusion principle applied to the system of subsets generated by k. For exam- 
ple: 7130 = F30—F10—F15+F5 and 77,60 = F^q—F^q— F20— F12+F10+F6+F4— F2. 



6 Disjoint sum of graphs 

The results of the previous sections may be formulated as follows. On the one 
hand we have a category Q of finite graphs. On the other we have the category 
of binary-generated trees T and the category of loop ensembles C: 

i 

The construction described in this paper is a map 

* : g — > C x T (6.2) 
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which associates to each graph G two elements, L{G) e C and T(G) £ T. The 
digraph of evolution G* is a semidirect product of these two, namely, if {L, °r>} 
and {T, S^} are viewed as digraphs 

G* = {L x T, ^} = L k T 

The problem is to find how various operations and maps in the category of 
graphs are reflected in the algebraic properties of the spaces C and T. 

From many such problems we answer a modest question concerning the evo- 
lution graph for a disjoint sum of two graphs, G\ U G 2 . Before we go on, a few 
remarks on the two categories: of trees and of loop ensembles. 

Universal family of binary generated trees. The family of trees T := gen {To ,h,h, ■ ■ ■} 

consists of all finite *-products of binomial trees. It forms a semigroup (T, *) 
satisfying: 

1) T 1 *T 2 =T 2 * Ti 

2) (Ti * T 2 ) * T 3 = Ti * (T 2 * T 3 ) (6.3) 

3) T * I =T 

Thus (T, *) forms a commutative free monoid with exponentials from N. 

Universal loop space. Similarly, let us introduce a "space" C over N spanned 
by the the set {L 1 ,L 2 , ■ . .}, where represents a loop length i. Using additive 
notation, a collection of loops (loop ensemble) will be denoted as a formal sum, 
e.g. L = L 3 + 2L 4 = {A, □, □}. Now we shall turn C into an algebra by intro- 
ducing a product * : C x C -> £, which for single loops L a and i b of length a 
and 6 respectively is defined 

L a *L b = (a\b) L [a \ b] (6.4) 

where we denote (a\b) = gcd (a, 6) and [a|fo] = 1cm (a, 6). 
Theorem 6.1 Tfe fn'p/e {£,+,*} satisfies 

1) L a * Lb = Lb * L a 1 ) L a + Lb = Lb + L a 

2) {L a *L b )*L c = L a *(L b *L c ) 2') (L a + L b ) + L c = L a + (L b + L c ) 

3) Li * L = L for any L 3') L + L = L for any L 

4) (L a + L b ) * L c = L a * L c + L b * L c 

Proof: Let us show item 2): Recall that (a\b) = Using (6.4) we get 

/r j.rAj.r a 'b T j. r a 'k [a|6]-c r o-6-c r 

^ * ^b) * -^c - pb] ^[a|6] * - pb] • pbH ^[o|6|c] ~ J^\b\c\ H^Wc] - 

(a\b\c) L[ |b|c]- Since [ | ] is an associative operation, property (2) holds. □ 
The definition in Eq. 6.4 is chosen to make the following look simple. 
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Proposition 6.2 The evolution digraph of a disjoint union of two graphs is 

(d + G 2 )* = {T(Gi) * T(G 2 ), L(Gi) * L(G 2 )} (6.5) 

jx>/zere f/ze ster operations are defined in (6.4) for loops, and (6.4) for the characteristic 
trees. 



7 Power Graphs 

Definition 7.1 The jfti power of a graph G is the graph G p determined by the jfi 1 
power of the harmonic matrix of graph G, a p . 

Clearly, G p has the same set of vertices as G, but a different structure of edges. 
That the powers of a graph are well-defined is obvious, since a p fl = a p ~ 1 (aft) = 
and (a p ) T = (a T ) p = a p . An interesting question is to determine the evolu- 
tion digraph of graph G p from the the evolution digraph of G alone. 

From the basic properties of the monoid M(a) we obtain the length of the tail 
and rank of the cycle group of M(a p ) to be 

tail: |7V)|=Ent(|T(a)|/g) 
loop: \L{a p )\ = \L(a)\ // q , 

where / denotes "funny division," defined for integers by a//b = [a\b]/b (recall 
that [a\b] = lcm (a, b)). For instance 6//1 = 6, 6/2 = 3, 6/3 = 2, 6/4 = 3, 
6/5 = 6, 6/6 = 1. 

Proposition 7.2 Let the evolution digraph G* of graph G have the following charac- 
teristic tree and the loop ensemble: 

T(G) =/N/N...*^ = fi 

i=l 

m 

L(G) = c\L\ + c 2 L 2 + . . . + c m L m = J2 CiLi , 



where k = \T(a) \ and m = \L(a)\. Then the evolution digraph (G q )* consists of 

(7.2) 



T(c) = n ir where w *= e 5- h 

i=l j=g(i-l)+l 



L(G") = E <k ■ (i\q) L l//q 



i=i 



where k' — \T(a p \ and m! = \T(a p )\ are those of (7.1) above. 

Proof: The proof of this is simpler in combinatorial terms. Express G* as the 
product of the loop ensemble and the characteristic tree. A state S is equiva- 
lent to a choice of a node in one of the loops in L(G), and a node in each binary 
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tree constituting T(G) (imagine them as "glowing points"). The action of the 
harmonic operator on S is represented now by a shift of each glowing point by 
one position: along the loop, and down each tree. Tracing their g-step move- 
ments brings the conclusion to the proof. ^ 

Remark: Every integer geN defines a linear map on the universal loop space 
q : C — > C which transforms each single loop L a into a sum of loops accord- 
ing to 

L a — > (a\q) L a//q 

as suggested by (7.2). For instance 4(L 6 ) = 2L 3 , which is due to the two distinct 
paths of "4-step jumps" in a hexagon. The second part of Proposition 7.2 may 
be now expressed simply 

L(G") = qL(G) 

Notice that operators a representation of natural numbers with multipliation, 
namely ab = a ■ b. On the other hand, we have the adjoint representation of the 

algebra {£, *}, where for a e N we define si t = L a * Li. Each of these repre- 
sentations is of course commutative, but the common algebra of both kinds of 
operators is not. For instance: 

f L 2 i if i = odd 

* A I A * 

2 2Li=< 2Li if 2|z but 4 /fi whereas 22 Li = 
{ AL l/2 if 4|i 

8 Conclusion 

We defined the harmonic operator of a graph as a topologically motivated en- 
domorphism of the space of states of graphs. The states (co-states) are meant 
as vectors in the linear space over Z formally spanned by the vertices (edges). 
The paper concentrates on the former. We defined a dynamical system on a 
graph as a recursive application of the harmonic operator. We found the types 
of evolutions (evolution digraphs) and classified them. An algorithm to fore- 
tell the evolution digraph for a particular graph has consists of Theorems 4.2 
and 5.2. 

We see it as a map 

g ^ CxT (8.1) 

which associates to each graph G two vectors, L{G) e C and T(G) e T. Trees 
and loops correspond to structure of invariant subspaces of a; any number de- 
rived from it may be viewed as an invariant of the graph. 

Map (8.1) is a morphism of objects, and the correspondence for operations of 
the disjoint union and for the graph power were shown. The problem of find- 
ing such a correspondence for other operations on graphs may prove to be 



f 2L t if i = odd 
\ 4L i/2 if i = even 
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rather difficult. 

The inverse problem asks whether a given digraph represents a harmonic evo- 
lution of some graph. Equivalently, the problem is to characterize the images 
L(G) C L and T(Q) C T. We offer these two conjectures. 

Conjecture 1: Any finite product of binary trees of height at least 1 is a char- 
acteristic tree for some graph. 

Conjecture 2: Denofe H = {2 l | i g N} = {1,2, 4, 8,...} (the set of "highly 
even numbers"). A formal sum L = is said to be admissible, L g La, if 

the following is satisfied: 

1) m g H 

2) i-rti g H for each p such that n p ^ 

i\p 

3) if Hi ^0 or rij 7^ 0, then ^ 

4) log 2 i • nj is even 

where we denote = l.c.m. (i, j). Then fhe sef of harmonic loop ensembles 
of finite graphs coincides with the above admissible set, C(Q) = Ca- 

Further exercises and problems: 1. rewrite the theory for the co-harmonic op- 
erator b = Sd acting on the edges and defining the harmonic "co-evolution" on 
graphs. Can one determine the corresponding co-evolution digraph from the 
evolution graph alone? 2. Replace the field Z 2 by another commutative ring 
(like Z 3 ). 3. Classify the rules for evolution on graphs other than harmonic 
rules. 
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